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of X-ray beam from rough surface 
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Abstract 

A new method for investigation of X-ray beam propagation in a rough 
narrow dielectrical waveguide is proposed on the basis of the numerical 
integration of the quazioptical equation. In calculations a model rough 
surface is used with the given distribution of hights of roughness and 
given correlation properties of the surface. The method is free from 
the limitations such as infinitcncss of the surface length and plane wave 
approximation which are nessesary for application of Andronov- Leon- 
tovich method. Our method can be applied to any surface with given 
nonhomogeniety and distribution of transitional layer. 
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In the middle of 1970-s the experiments on chanalization of soft X- 
rays |Q and middle range X-rays |2| in glass tubes and also for filtering 
of hard part of X-ray radiation in bent tubes bundles [|| . The workabil- 
ity of such systems was demonstrated for transport of X-ray radiation. 
Later the interest to this manner of X-ray steering was connected with 
the suggestion to use the samples of many specially bent tubes for geo- 
metrical focusing and concentrating of X-ray beam Q| . The experiments 
were performed for transmission of X-ray beam through " gap-less colli- 
mator" - micron gap betwean two tightly pressed together glass plates 
[|[ ^| . The revival of interest to localizing of X-rays in thin glass capil- 
laries of changing diameter in the condition of total external reflection 
at grazing incidence on smooth surface was connected with the attempts 
of microfocusing in a narrowing tube [Q (till 7/j,m) or in a narrowing 
polycapillary system consisting of a great number of melted- in to- 
gether thin glass tubes forming one hexahedral block. The problems of 
the structure of the focus of such system deals with the work |j| . 

The channalisation of X-ray radiation in hollow glass tubes is pos- 
sible due to the effect of total external reflection ]ic|j . The dielectric 
permeability of substance for electro-magnetic radiation with the energy 
higher than binding energy of electrons in atoms can be approximately 
written with the use of plasma frequency uj v 

e(uj) = 1 - (lu p /uj) 2 . 

Because in X-ray energy range e(a->) < 1 and optical density of glass 
is smaller than the vacuum (or air) optical density the effect of total 
external reflection takes place when the rays come from outside and for 
this reason is called total external reflection (TER). 

When X-ray radiation comes on sufficiently smooth surface in the 
conditions of TER the radiation penetrates only at small depth of the 
order of 60 A, that maintains the effective reflection and channalisation 
of radiation in hollow glass tubes. 

Besides of the absorption the effectiveness of reflection depends strongly 
on the scattering of radiation on roughness of the surface. The most of 
the results on the influence of roughness on the reflection of X-ray ra- 
diation from surfaces were obtained in the approximation of Andronov 
and Leontovich |TT| . The revue of the results can be found in [ 
approximation is based on the suggestion that the initial wave coming on 
the surface is a plane wave and the scattered wave can be found for small 
perturbations of the boundary between two media. That is why the in- 
vestigation of scattreing of wave coming on absorbing surface at small 
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angles is interesting because in this case the Andronov and Leontovich 
approximation apparentely can not be applied any longer. 



1 The derivation of master equation. 

Maxwell equations are |D| 

_ _ edE 471-cr 
rotH= -— + E, 

c at c 

rotE = -^ — , 

c at 

divE = 0, 
divH = 0. 

We can exclude the vector H by applying to the second equation 
the operation rot 

rotrotE = graddivE - AE. 
With account of the third Maxwell equation we get wave equation: 

e/i<9 2 E 47r<7 ii (9E 
Furthermore, represent the vector of the electric field in the form 

E=^A(x,y,z)exp[i(ut-kz)] + .. (2) 

k = v^o^- P 11 tlh s case eo is dielectrical permeability of air, e- dielec- 
trical permeability of glass.) It can be shown that the speed of changing 
of the amplitude A(z, x, y) in the beam is different for the transversal 
coordinates x,y and longitudinal coordinate z ]l3}] . The evolution along 
the direction z is much slower than along transversal coordinates. Thus 
after substituting of (2) into wave equation (1) we can neglect terms 
d 2 A/dz 2 in comparison with k(dA/dz) and d 2 A/dx 2 , resulting to so 
called parabolic equation of quazioptics: 

2ik— = Aj_A + k 2 — A (3) 
oz eo 

where the difference between the dielectrical permeabilities Se = e — e% 

depends on coordinates and includes imaginary part corresponding to 

absorption in the substance. 
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2 The method of simulation of rough sur- 
face. 



Under the total external reflection at grazing incidence of X-ray beam 
the most important is scattering within the plane of incidence because 
the scattering withing the plane of the interface is small Jl2|. First of 
all it is due to the fact that the wavelength of the radiation is much 
smaller than the characteristic scale of inhomogeneities. That is why 
for simulation of scattering of X-rays at TER conditions it is enough to 
account for only scattering in the plane of incidence within 2-dimensional 
model. In this case the value Se in the right side of equation (3) becomes 
the function of coordinates z and x : 



Se(x, z) 



0, x > £(z) 
e - e , x < £(z) \ ' 

where the function £ (z) is the hight of the rough surface profile and can 
be regarded as a random value (see Fig.l()). 

Stationary random value on the interval (0, Z) can be expanded in 
the Furier specious 

oo 

£(z) = Vk cosu>kZ + Uk siau>kZ, < z < Z, 
fe=i 

where Vk Uk are random amplitudes of the harmonics LUk = ku>x, 

u>i = 2n/Zi 

Z\ is maximum spacious period of random sequence. In the discrete 
representation 

kirn rr . kim 



= 2^VkCos—— + U k sin—-, n = l,N; 

k=0 

where Vk and Uk uncorrelated random numbers with the dispersion u\ 
and zero mean value. The dispersion of the harmonics is 
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3C: 



i?(0 cos(fcwiOde, (4) 



o 



where is correlation function of the random value £(z). For normal 
random processes the amlitudcs Vk and Uk must have normal distribu- 
tion G3]. The expression for ^(n) can be also repesented as 



3 



V—v Kirn . 



(5) 



where Ek is random coefficient with Rayleigh distribution with the pa- 
rameter a equal to er^ , where otk is random phase of the harmonics with 
the uniform distribution on the interval (0, 2ir). 

For choosing of the number of harmonics m a relationship can be 
used 



so the summ of dispersions <j\ must be equal to the dispersion of the 
simulated process. 

3 Numerical results and discussions. 

The numerical method for solving of the equation (1) was used earlier for 
study of motion of channeled electrons in single crystals |l5| and reflec- 
tion of positrons from slanting cut single crystals pl| . The method im- 
plies step by step calculation of the amplitude of the X-ray wave A(x, z) 
along the surface of reflection starting from its initial value at z=0. 

For simulation of a rough surface the code was constructed calculat- 
ing random sequence of numbers in accordence with representation 
(5). With the given correlation function by relationship (4) the ampli- 
tudes Ek of spectral components were determined. For simulation of 
random phase otk random numbers generators were used. Fragments 
of random surface £(n) is presented on on Fig. 1(a). With account of 
several realisation of numerical process with different number of points 
correletion functions were calculated again along with propability den- 
sity distribution. On Fig. 1(b) the comparison of given correlation 
function (curve 1) and correlation functions calculated with realisation 
of n=1000 points (curve 2) and n=10 000 points (curve 3 that nearly 
coinsides with 1). It is obvious that with the increase of n correlation 
function of the process £(n) approaches the given function. It was also 
shown that the distribution of the propability density for random pro- 
cess £(n) approaches Gaussian distribution when the number of points 
in the realisation increases (Fig. 1(c)). 
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So we can affirm that the created code simulates random surface with 
the given statistical characteristics and Gaussian propability density dis- 
tribution. 

The calculations were made for the radiation with the energy 10 keV 
and for interaction with rough glass surface. The width of the correlation 
function was chosen 5/im. The angle of total external reflection is tp c = 
3 • lQ- 3 rad. 

When the angle of incidence of X-ray wave is not zero the surface can 
be regarded as infinite and incident wave as a plane wave. On Fig. 2 the 
distribution of the intensity of the radiation along the coordinate x (the 
axis x is perpendicular to the surface) for reflection from smooth surface 
for the angle of incidence ip = 10~ 3 rad. The range < x < 2100A 
corresponds to the substance layer. The oscillations of the intensity 
are caused by the interference of incident and mirror reflected waves. 
At some distance from the surface that is defined by numerical scheme 
parameters the amplitude was smoothly truncated that is dependent on 
the use of the Furier transformation on the x coordinate for solving the 
equation (3) (see Fig. 2). 

With the increasing of roughness amplitude the angular spectrum of 
the reflected beam the mirror reflected beam decreases and at the same 
time numerous random maxima arise so that clear interference diffrac- 
tion picture disappears near the surface. Simple estimate for roughness 
tolerance of total external reflection observation was given in |l7j 

h < Xo/Sipo, 

where Ao is the radiation wavelength. As a clear illustration of inverse 
dependence of tolerable roughness amplitude on the incidence angle can 
be used angular spectra of reflection shown on Fig. 3. The spectra were 
calculated for angles 0.5 • 10~ 3 rad. (a),10~ 3 rac?.(b), 2 • 10~ 3 rad.(c) for 
the same roughness amplitudes 200A (here and further the mean squared 
roughness amplitude is used) . In the first case the influence of roughness 
on the spectrum is small but in the last case the mirror reflected peak 
can not be observed in practice. 

The calculations of the distribution of the intensity of radiation near 
the smooth reflecting surface can show the depth of penetration of X- 
rays into the surface under the conditions of TER as a function of the 
incidence angle and for <p = 10~ 3 rad. it is 60A, for tp = I0~ 4 rad. - 10 A. 
For reflection from rough surfaces the depth of penetration (counted 
from the mean value of the roughness) is approximately equal to the 
mean squared amplitude of roughness. 



5 



With the created complex of codes calculations were performed for 
transmission of X-ray radiation in a 2-dimcntional rough capillary for 
zero and near to zero entrance angles. Transversal dimensions of cap- 
illaries were chosen 1 — 2/ira, that caused diffractional spreading of the 
plane wave to the values 0.5 • 1CP 4 — I0~ 4 rad. The length was chosen as 
1 — 2cm. 

On Fig. 4 the dependence of integral intensity (over the transversal 
dimension of the capillary) on the distance from the entrance to the 
capillary. The angle of incidence ipo = 0, the width of the capillary 
is 2/um, the roughness is 1200A. Abrupt falling of the intensity of the 
radiation in the input of the capillary is connected with the absorption 
of the radiation coming on glass butt-ends and transmitting within the 
substance. This part of the radiation can be regarded as fully absorped 
at the distance z < lOOO^m from the input of the capillary. The analysis 
of the distribution of losses along the capillary length for z > 1000/im 
reveals its considerable dependence on the distance from the capillary 
input. The deeper radiation penetrates the capillary the smaller loss 
normalized to the length unit is that is apparently related to consequent 
decay of modes the most strongly penetrating into the substance. 

The losses of radiation in the capillary are the greater the the greater 
capillary wall roughness is and the bigger the angle of incidence of the 
wave into the capillary is. The partial dependence of losses of the radia- 
tion (size for zero incidence angle and angle value 1CP 4 rad is shown on 
Fig. 5. The losses were calculated on the range 0.25cm < z < 1cm from 
the capillary input. It is worth noting that losses are not zero even for 
the capillary with smooth walls. 

The dependence of losses of radiation in capillaries of various widths 
with wall roughness 800A and 400A and also without roughness for zero 
incidence angle are shown on Fig. 6. The increase of losses when the width 
of the capillary Ax is decreased can be accounted for by the broadening 
of the plane wave incident on the input of the capillary A(p—Xo/ Ax. 

On Fig. 7 radiation angular spectra are shown for capillaries with 
the length 1 cm of various widths with smooth and rough walls. The 
calculations were made for zero incidence angle. The decrease of the 
square under the spectral peak along with the increase of roughness can 
be accounted for by higher absorption of the radiation but the shape 
and the width of the spectral peak practically does not depend on the 
roughness size. 
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4 Conclusion. 



So the created code can simulate interaction of X-rays with rough surface 
(including zero gliding angle) without supposition of infinite plane wave 
but with direct account of given surface relief. That is the advantage 
of the approach over analytical methods (incliding Andronov-Leontovich 
approimation). In comparison with numerical methods using geometrical 
optics the proposed method takes into account the wave nature of the 
radiation (diffraction) . 

The obtained results enable also to give new interpretation of the 
experiments with "gapless collimator" published in 1981-1984 J5|, |). Ac- 
coding to the data of angular distanse between interference peaks in 
transmitted beam with account of diffraction the gap width is 10/im but 
not 1 fiui as it follows from || without account for diffraction. In this 
case different number of transversal modes is excited depending on the 
tilt angle of input beam. On the other hand if in glass plates unpolished 
transversal band is left the independence of output beam on the input 
beam tilt can be explained due to strong absorption of higher transversal 
modes along with retaining of lowest symmetric mode (see Fig. 7). 
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Fig. la. Fragment of calculated surface with account of equation (5). 
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Fig. lb. 1 - given initially correlation function; 

2(dashcd) - correlation function calculated with realization of random se- 
quence (n) with 1000 points, 

3 - correlation function for sequence with 10000 points. 
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Fig.lc. 1 - gaussian distribution with given parameters, 

2 - distribution of density probability calculated calculated with realization 
of random sequence (n) with 10000 points. 
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Fig. 2. Distribution of intensity of X-ray radiation in the vicinity of smooth 
reflective surface in the conditions of the Total External Reflection at the 
angle of incidence on the surface ipo = 10~ 3 rad. Axis x is perpendicular 
to the surface. Reflecting layer of the substance is within < x < 2100A; 
reflecting surface is marked with dash line. 
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Fig. 3. Angular spectrum of reflection of X-ray radiation from the surface with 
roughness amplitude 200A for incidence angle to the surface <po = 0.5 • 10~ 3 
rad. 
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Fig.3bc. The same as Fig. 3 for ipo = 10 3 rad (b) and tpo = 2 • 10 A rad 
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Fig.4. Dependence of integral over the width of the capillary x = 2 fim 
intensity of radiation on the depth of capillary Z. Angle of incidence of plane 
wave to the input of the capillary tpo = 0, wall roughness is 120CL4. 
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Fig.5. Dependence of losses of radiation (%) calculated at the distance 
0.25cm < Z < 1.0cm from the input of the capillary on the averaged am- 
plitude of roughness at the capillary walls (a 2 ) 1 / 2 - mean squared deviation 
for roughness {;■ Angles of incidence of X-ray beams to the input of capillary 
<po = 10~ 4 rad(l) ipo = rad (2). The capillary width is 2/im. 
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Fig. 6. Dependence of losses of radiation (%) calculated within the range 
0.25 < z < 1.0cm from the input of the capillary. Angle of incidence ipo = 0. 
Roughness of capillaries is 80()A(1), 400A(2), smooth walls (3). 
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Fig. 7a. Angular spectra of the radiation at the output of capillaries with 
length 1cm for smooth capillaries (solid lines) and with rough capillaries 
400A(dashed lines). Width of capillaries is 0.5/wra. Angle of incidence is 
<P0 = 0. 
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Fig. 7b. The same as Fig.7a but for roughness 800A and width of capillaries 
1/im. 
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